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In this paper, we calculate the diffusion coefficients that are related to the neutrino opacities con- 
sidering the formation of nuclear pasta and homogeneous matter at low densities. Our results show 
that the mean free paths are significantly altered by the presence of nuclear pasta in stellar matter 
when compared with the results obtained with homogeneous matter. These differences in neutrino 
opacities certainly influence the Kelvin-Helmholtz phase of protoneutron stars and consequently the 
results of supernova explosion simulations. 
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I. INTRODUCTION 

When massive stars (8 Mq < M < 30 Mq ) exhausts 
its fuel supply, the forces that support the stars core 
quickly retreat, and the core is almost instantly crushed 
by gravity, which triggers a type II supernova explosion. 
The remnant of the gravitational collapse of the core of a 
massive star is a compact star or a black hole, depending 
on the initial condition of the collapse. Newly-born pro- 
toneutron stars (PNS) are hot and rich in leptons, mostly 
e~ and v e and have masses of the order of 1 - 2 M 
During the very beginning of the evolution, most of the 
binding energy, of the order of 10 53 ergs is radiated away 
by the neutrinos. 

The composition of protoneutron and neutron stars 
remains a source of intense speculation in the litera- 
ture. Whether their internal structure is formed by nu- 
cleons and leptons, by other light baryons and leptons, 
by baryons, leptons and quarks (bearing or not a mixed 
phase), by baryons, leptons and kaons or by other pos- 
sible composition, is still unknown. The neutrino-signals 
detected by astronomers can be used as a constraint to 
infer protoneutron star composition [2, 3|. For the same 
purpose, theoretical studies involving different possible 
equations of state obtained for all sorts of matter com- 
position have to be done because the temporal evolution 
of the PNS in the so-called Kelvin-Helmholtz epoch, du- 
ring which the remnant compact object changes from a 
hot and lepton-rich PNS to a cold and deleptonized neu- 
tron star depends on two key ingredients: the equation 
of state (EoS) and its associated neutrino opacity at su- 
pranuclear densities 0, HI • 

Neutrinos already present or generated in the PNS hot 
matter escape by diffusion (not free streaming) because of 
the very high densities and temperatures involved. The 
neutrino opacity is calculated from the scattering and 
absorption reactions that take place in the medium and 
hence, related to its mean free path, which is of the or- 
der of 10 cm and much smaller than the protoneutron 
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star radius 5]. In the diffusion approximation used to 
obtain the temporal evolution of the PNS in the Kelvin- 
Helmholtz phase, the total neutrino mean free path de- 
pends on the calculation of diffusion coefficients, which, 
in turn, depend on the chosen EoS. At zero temperature 
no trapped neutrinos are left in the star because their 
mean free path would be larger than the compact star 
radius. 

A complete equation of state capable of describing 
matter ranging from very low densities to few times sa- 
turation density and from zero temperature to around 
50 MeV is a necessary step towards the understanding of 
PNS evolution. The constitution of the PNS crust plays 
a definite role in the emission of neutrinos. For this rea- 
son, the pasta phase, present in very low nuclear matter 
as the crust of PNS are included in the investigation of 
the neutrino opacity in the present work. 

A few words on the pasta phase follow. It is a frustra- 
ted system |gjlO| present at densities of the order of 0.006 
- 0.1 fm -3 [111 in neutral nuclear matter and 0.029 - 0.065 
fm -3 [l2l . Il3l ] in /3-equilibrium stellar matter, where a 
competition between the strong and the electromagnetic 
interactions takes place. The basic shapes of these com- 
plex structures were named |6J after well known types of 
cheese and pasta: droplets = meat balls (bubbles = Swiss 
cheese), rods = spaghetti (tubes = penne) and slabs = 
lasagna, for three, two and one dimensions respectively. 
The pasta phase is the ground state configuration if its 
free energy per particle is lower than the corresponding 
to the homogeneous phase at the same density. 

The evolution of PNS and simulation of supernova ex- 
plosion have already been considered for different matter 
compositions, some with the inclusion of the pasta phase 
d Ht 0, 03, El]. From 0, 13 one can see that the trans- 
port properties are significantly affected by the presence 
or absence of hyperons and of the mixed phase in hy- 
brid stars. In [7| the linear response of the nuclear pasta 
to neutrinos was calculated with a semi-classical simu- 
lation and the muon and taon neutrinos mean-free path 
were described by the static structure factor of the pasta 
evaluated with Metropolis Monte Carlo simulations. In 
14] rod-like (two dimensions) and slab-like (one dimen- 
sion) pasta structures were included in the calculation 
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of neutrino opacity within quantum molecular dynamics. 
A very interesting conclusion was that the pasta phase 
occupies 10-20% of the mass of the supernova core in the 
later stage of the collapse. 

In the present work we investigate the influence of the 
pasta phase on the neutrino opacity by showing the ef- 
fects on the diffusion coefficients. The pasta phase is 
calculated with the coexistence phases method (CP) in a 
mean field approximation [IJ [lj, EU • We consider only 
nucleons and leptons in the EoS in /3-equilibrium. In the 
pasta structure only electron neutrinos are considered. 

In section II we present the formalism used to obtain 
the equation of state, in section III the recipe used for the 
construction of the pasta phase is outlined, in section IV 
the expressions used to calculate the neutrino cross sec- 
tions and related mean free path are given and in section 
V our results are shown and the main conclusions are 
discussed. 



II. FORMALISM 

We consider a system of protons and neutrons with 
mass M interacting with and through an isoscalar-scalar 
field <j> with mass m s , an isoscalar- vector field with 
mass m v and an isovector- vector field h p with mass m p 
described by the well known non-linear Walecka model 
(NLWM) [lTj ■ We impose /3-equilibrium and charge neu- 
trality with neutrino trapping at finite temperature. At 
zero temperature no neutrinos are left in the system. 

The Lagrangian density reads 

C= A'+Ar + £ U +£ P + A, (l) 

j=p,n l=e,v 1 



where the nucleon Lagrangian reads 

£ i= ^.[7^-M*]^, (2) 
were M* = M — g s (f> is the effective baryon mass and 

W = id* - g v V - ^-t • b". (3) 

The meson Lagrangian densities are given by 

L* = \ ta W - ™l<? - \*4>* - ^<t>*) , (4) 
C^J-ifl^lmXp), (5) 
C p = ±(-±B^-B^ + rn*h ll -V>y (6) 

where fi M „ = <9 M V„ - d u V p and B M „ = <9 M b„ - 9„b M - 
g p (ba x by). The lepton Lagrangian densities read 

A =i>i h^id" -mi]i>u (7) 

where m e is the electron mass and the neutrino mass is 
m„ = 0. 

The parameters of the model are three coupling cons- 
tants g s , g v and g p of the mesons to the nucleons, the 
nucleon mass M, the electron mass m e , the masses of 
the mesons m s , m v and m p and self- interacting coupling 
constants n and A. The numerical values of the parame- 
ters used in this work and usually referred to as NL3 [lH 
are shown in table |U They are fixed in such a way that 
the main nuclear matter bulk properties are the binding 
energy equal to 16.3 MeV at the saturation density 0.148 
fm -3 , the compressibility is 272 MeV and the effective 
mass at the saturation density is 0.6 M. 



Model 


9s 


9v 9 P 


M m e m s m v m p k/M A 


NL3 


10.217 


12.868 8.948 


939.0 0.511 508.194 782.501 763.0 4.377 -173.31 



Tabela I. Parameters set used in this work. All masses are 
given in MeV. 



From the Euler-Lagrange formalism we obtain the meson fields by their mean values 
equations of motion for the nucleons and for the meson 
fields: 

V72^ 2,1 ,2, 1 iA 3 /q\ <l> (4>) = 4>0, (11) 

V* = 2* +3!^ ~9sPs, (8) V^(V,) = V , (12) 

V 2 V = mlVo - g vPBl (9) 6 M _). (^) = b , (13) 
V 2 & = m>o - %Pz, (10) 



where p s , pb and p% are defined next. By replacing the the equations of motion read 
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Vn = 



2mV° 6mV° 



tPb, 



bo = TT^PS' 

2m~ 



(14) 
(15) 
(16) 



where ps — p P + p n is the baryonic density and p3 
Pp — Pm Pp and p n are the proton and neutron densit 
given by 



ics 



Pi = 2 J (^rT3^' + ~ 3 =P> n (17) 



where f j± = 1/(1 + exp [(e, T Vj)/T\), ej = Vp 2 + M* 2 
and Vj = Hj — g v Vo — g p T^bQ, where T3 is the appropriate 
isospin projector for the baryon charge states and pij are 
the nucleon chemical potentials. The scalar density p s is 
given by 



^ f d s p M* 



+ fj-)- a* 



where the electron and electron neutrino energy densities 
are 

£l = 2^j dp P 2 \fp 2 + m Kfi+ + fl~)> ( 28 ) 
and electron and electron neutrino pressure are 

P l = JL [dp -^=(f l+ + /,_). (29) 



6tt 2 



Vp 2 + ' 



The electron density p e and electron neutrino density p v 
are given by 



(30) 



where g e = 2, g v = 1, fi± = 1/(1 + exp [(ej T Pi)/T]) 
with ei = yj p 2 + mf and \x e is the electron chemical po- 
tential, e„ is the electron neutrino energy and p, v is the 
electron neutrino chemical potential. The condition of /3 
equilibrium in a system of protons, neutrons, electrons 
and trapped electron neutrinos is 



Pp= Pn~ Pe+ Pv 



(31) 



We impose neutrality of charge as p p = p e and fix the 
lepton fraction 



The thermodynamic quantities of interest are given in 
terms of the meson fields. They are the total energy 
density 



(19) 



l—e,u 



with 



f = ^ E / d P P 2 Vp 2 + M* 2 (f 3+ + /,•_) (20) 



^v 2 +^b 2 +4ti 



6^24^' 



the total pressure is 



Pt = P+ E P << 



with 



p =i.S /+■ 



l—e.u 



,2 p 



yjp 2 + M* 



2 Vo+ 2 b ° 2 00 



(21) 
(22) 

(23) 



f (/ J+ + /i_) (24) 

(25) 
(26) 



and the total entropy density 

5 = i(fir + - HPi - E MiPi). (27) 



Yr. 



PS 



(32) 



Notice that muons are not considered in the present 
calculation. 



III. COEXISTING PHASES: NEUTRAL 
NUCLEAR MATTER WITH NEUTRINO 
TRAPPING 



The formation of pasta phase has been studied lately 
with great interest 0, [l9| . Next we show the main steps 



for the calculation of the pasta pha se with the coexistence 
phases method based on [20I |2l|. For further details, 
please refer to [ill G3 ■ 

For a given total density ps and lepton fraction Yl we 
build pasta structures with different geometrical forms in 
a background nucleon gas with f3 stability and neutrino 
trapping. This is achieved by calculating from the Gibbs 
conditions the density and the particle fractions of the 
pasta and of the background gas so that in the whole we 
had to solve simultaneously the following seven equations 

P\vl vl M* 1 ) = P"(!tf, itf, M*Y33) 
1 - ■■" (34) 
(35) 
(36) 



Pv Pv J 



Pn 
Pi 



2 1 T t j 

m s<Po + 77 OA 



I\2 



2 J,IJ 



m s (j) 



II\2 



^ ■ II\3 



9sP s 



l—e.v 



f(p I v -pl) + (l-f)(pl I -pl I )=^ 



(37) 

(38) 
(39) 
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where I and 77 label each of the phases, / is the volume 
fraction of phase / 



P 1 



ii 



ii 



(40) 



The total pressure is given by Pr — P 1 + P e + Pv ■ The 
total energy density of the system is given by 



£ = f£ T + (1 - f)£" + £ e + £„ + £ surf + £ C ouu 
with £ sur f = 2£ Coul [H, [H] , and 
2a 



(41) 



£, 



Coul 



42/3 



{e^Y'^aDipl-pl 1 )?/*, (42) 



where a = f for droplets, rods and slabs and a = 1 — / 
for bubbles and tubes, a is the surface energy coefficient, 
D is the dimension of the system. For droplets, rods and 
slabs, $ is given by 



D-2 



/-l-ln(/) 
D+2 • 



+ / 



1 

D+2 



D = 1,3, 



(43) 



£ = 2, 



and for bubbles the above expressions are valid with / 
replaced by 1 — /. The surface tension plays a signifi- 
cant role on the appearance of the pasta phase. In our 
treatment of the surface tension we essentially follow the 
prescription given in [ill . [l3|, but some comments on 
the importance of the surface energy on the calculation 
of the pasta phase are mandatory. It has been shown that 
the existence of the pasta phase as the lowest free energy 
matter and of its internal structures essentially depends 
on the value of the surface tension [t| [ll], [l3|, [l(| [2J| . In 
the present paper the surface energy coefficient is para- 
metrized in terms of the proton fraction according to the 
functional proposed in [25(, obtained by fitting Thomas- 
Fermi and Hartree-Fock numerical values with a Skyrme 
force. The same prescription was used in [Til. [l3l] . Howe- 
ver, a better recipe is to consider the surface energy coef- 
ficient in a consistent way, in terms of relativistic models. 
In the surface energy was parametrized according to 
the Thomas-Fermi calculations for three parametrizati- 
ons of the relativistic NLWM. The Gibbs prescription 
was used to obtain the a coefficient which is the appro- 
priate surface tension coefficient to be used [2(| [27j ■ This 
improvement will be added to our calculations in a forth- 
coming work. 



IV. NEUTRINO CROSS SECTIONS 

To calculate neutrino opacities and mean free paths we 
consider [f| neutral current scattering reactions 



■ n — t v e + n, 
\-p -)• u e + p, 



(44) 
(45) 



and charged current absorption reactions 



n 
P 



- n. 



The cross sections for reactions J44 
(|T7l) employed in this study follow [5] : 
Reaction ([ti]): 



(46) 
(47) 

461) and 



On = < 



/ \ 2 

<J ref — (x) ' 

(*) (^) , 

°~s K) (^ ' 

i (£) (&) m 



nND, vD or uND, 
nD, vND, [28], 



nD, uD, [29, 30]. 



Reaction (H5 



pND, vD or i/ND, 



M , pD, vND, 
o-„ ( , pD, vD, [29] 

Reactions (HB]) and (l47l): 



(48) 



(49) 



On = < 



a ref (l + 3g%), nND, vND, 

a ref (l + ig\) , nND, vD or vND [31], 

o ref (l+3g 2 A )(l) (^)x 
0, nD, Y L < 0.08. 



nD, vD, [281 



(50) 

In this expressions p f and e f mean the Fermi momen- 
tum and Fermi energy of the degenerate neutron. Y e , 
Y n , Y p , Yl, are the electron, neutron, proton and lepton 
fractions. ND denotes the non degenerate regime, while 
D denotes the validity in case of degenerate particles. 
(T = 1-76 x 10~ 44 cm 2 and qa — 1.254. Regions of in- 
termediate degeneracy are also handled: the degenerate 
and non degenerate sectors for both the baryons and the 
neutrinos of the cross sections detailed in equations (|48[) , 
(|4"9"| and (fST))) are joined by a simple interpolation algo- 
rithm as done in [2|. The total neutrino mean free path 
in dense matter is given by 



K = 



1 



PnO r , 



PpOp + PBO a 



(51) 



Rosseland neutrino mean free paths are related with dif- 
fusion coefficients Dj [J] by 



A 



D k 



J™de„ e£/„(l - f v ) ' 



where 



POO 

D k = / de v e k v X v f v {l - f v ). fc = 2,3,4 
Jo 



(52) 



(53) 
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Figura 1. free energy per particle with the NL3 parametriza- 
tion obtained for T = 5MeV and Y L = 0.4. 




P (fm- 3 ) 



All contributions from neutrino opacities are related with 
the diffusion coefficients and can to be used as input 
to the solution of the transport equations in the equi- 
librium diffusion approximation to simulate the Kelvin- 
Helmholtz phase of the protoneutron stars [HJ ■ 



V. RESULTS AND CONCLUSIONS 

Before we tackle the problem of the consequences of 
the pasta phase on the diffusion coefficients, we display 
a characteristic figure of the free energy for the homoge- 
neous and pastalike matter obtained for T = 5 MeV and 
Yl = 0.4 in Fig. [TJ A similar figure is presented in Fig. 
9 of [16j . but obtained with a relativistic surface energy. 
One can see that the pasta phase ends when the free 
energy density reaches the curve for the homogeneous 
matter. Actually, at this temperature, the pasta phase 
interpolates between two regions of homogeneous matter, 
which is the preferential ground state at extremely low 
densities, as seen in Fig|T] Moreover, the size of the pasta 
phase decreases with the increase of the temperature and 
eventually, it no longer exists. It is also worth mentio- 
ning that neutrino free matter in /3-equilibrium presents 
a pasta phase smaller than matter with trapped neutri- 
nos [l3|, HH as a consequence of the fact that the latter 
presents a larger fraction of protons. According to stu- 
dies on binodals and spinodals underlying the conditions 
for phase coexistence and phase transitions (111 l33l . |34| , 
non-homogeneous matter with trapped neutrinos is ex- 
pected to be found until temperatures around T = 12 
MeV, depending on the model considered. 

We next show the diffusion coefficients Di , D3 and D4 
as function of the baryon density for different temperatu- 
res obtained for both homogeneous matter and the pasta 
phase. According to [TT|, El the densities where mat- 
ter becomes homogeneous depend on the proton fraction 



and on the temperatures involved, but it is always smal- 
ler than 0.1 fm" 3 for the NL3 parametrization and for 
the a values we consider in the present work. 

In obtaining the diffusion coefficients, the EoS was 
calculated as a grid where temperature ranges are in 
between and 50 MeV and densities vary from 0.005 to 
0.5 fm~ 3 . In our codes we have implemented a prescrip- 
tion given in [35| to evaluate the Fermi integrals so that 
the same codes run from zero (10 -9 ) to high temperatu- 
res. We have calculated the diffusion coefficicntcs only 
for baryonic densities above 0.005 fm~ 3 because the in- 
tegrals of type (|5"3l are very difficult to converge at lower 
densities. We show results for lepton fractions equal to 
0.2 and 0.4 because those are typical values necessary in 
the numerical simulation of protoneutron star evolution. 

In all figures the diffusion coefficients obtained with 
homogeneous matter join the curves obtained with the 
pasta phase at densities higher than the ones shown. For 
D2 calculated at T=5 MeV and Yl — 0.4, for instance, 
they cross each other at p = 0.12 fm -3 . Our codes inter- 
rupt the calculation once homogeneous matter becomes 
the ground state configuration, as depicted in Fig. [1] 
This means that there will always be a gap in the dif- 
fusion coefficients when the transport equations are cal- 
culated with the inclusion of the pasta phase. The same 
behavior is found at the pressure values for homogeneous 
and pasta phases at the transition density. 

From figures [2j [3] and |4] we can see that only three 
structures are found inside the pasta phase for the present 
model: droplets, rods and slabs as far as Yl = 0.4. For 
Yl = 0.2 only the first two structures remain. While the 
diffusion coefficients obtained with homogeneous matter 
is always smooth and continuous, a common trend of all 
the diffusion coefficients obtained with the pasta phase 
is a kink at very low densities in between 0.01 and 0.015 
fm~ 3 . The interpolation procedure we use depends on 
the quantities rji = (//$ — M*)/T,i = p, n. Whenever 
either r\ v or r\ n inverts its sign, these kinks appear, i.e., 
they are the result of the effective nucleon mass being 
greater than the corresponding chemical potential. Mo- 
reover, the pasta phase diffusion coefficients are always 
lower than the corresponding coefficients obtained with 
homogeneous matter. 

Our results for the diffusion coefficients D 2 and D4 are 
one order of magnitude larger than the results obtained 
in [32| . This difference can be explained because in the 
present paper all diffusion coefficients are calculated at 
very low baryonic densities. For larger densities the re- 
sults coincide. 

In summary we point out that in the present paper we 
have investigated the influence of the pasta phase on the 
neutrino opacity by calculating the diffusion coefficients. 
The homogeneous EoS was obtained with the NL3 para- 
metrization of the NLWM in a mean field approximation. 
The pasta phase was obtained with the coexistence pha- 
ses method (CP). 

Recent calculations for the pasta phase within the 
Thomas-Fermi approximation at finite temperatures [36j 
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show that the internal pasta structure is much richer as 
compared with the CP method we have employed in the 
present work. Hence, the dependence on the structure 
of the pasta phase is also of interest and this calculation 
is planned for different parametrizations of the NLWM. 
More sophisticated matter, which includes the a-particlcs 
should also be considered flU ]. 

We have checked that the neutrino interactions in 
warm and low baryonic densities with pasta formation 
show significant differences when compared with homo- 
geneous matter. Next the temporal evolution of the PNS 
will be calculated and, in face of the present results, we 
expect that the cooling and deleptonization eras will be 



influenced by the presence of the pasta phase at low den- 
sities. 

An obvious improvement is the inclusion of hyperons in 
the EoS. However, the pasta phase can still be computed 
just with protons, neutrons and light clusters because 
hyperons are expected to appear only at densities where 
the pasta phase is no longer present. 
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Figura 2. Diffusion coefficient D2 as function of baryon den- 
sity for different temperature and proton fraction values for 
homogeneous matter and pasta phase. 
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Figura 3. Diffusion coefficient D3 as function of baryon den- 
sity for different temperature and proton fraction values for 
homogeneous matter and pasta phase. 



14 
12 
10 
8 
6 
4 
2 





0.01 



0.02 



Rods 
Droplets 

Homogeneous matter 



0.03 



0.04 0.05 

P (fm" 3 ) 





14 


E 






12 - 


> 








5 


10- 






O 


8- 




CO 




Q 


6- 




4- 




2 




0- 




_ . Slabs 

Yl=0.4 — Rods 

Droplets 

NL3 Homogeneous mat te r| 



0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 

P (fm" 3 ) 




9 



Figura 4. Diffusion coefficient D4 as function of baryon den- 
sity for different temperature and proton fraction values for 
homogeneous matter and pasta phase. 





